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Within the weak-coupling BCS scheme we derive a general form of the coefficients in 
the Ginzburg-Landau expansion of the free energy of a superconductor for the case of a 
Fermi level close to a van Hove singularity (VHS). A simple expression for the influence 
of the VHS on the specific heat jump is then obtained for the case where gaps for different 
bands are distinct but nearly constant at the corresponding sheets of the Fermi surface. 
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1. Introduction 

The influence of a van Hove singularity (VHS) on the properties of supercon- 
ductors is a largely debated problem in physics of supcrconductivity.EJ Thus, we 
shall attempt here to analyze the influence of VHS in the density of states (DOS) 
on the jump of the specific heat AC at the critical temperature T^. This problem 
reduces to finding the coefficients of the Ginzburg-Landau (GL) expansion of the 
free energy. Particular attention will be paid to the account of multigap effects. 
Let us recall that 44 years ago Moskalenkoi predicted the existence of multigap 
superconductivity, in which a disparity of the pairing interaction in different bands, 
such as the s and d bands in transition metals, leads to different order parameters 
and to an enhancement of the critical temperature. Subsequently multiband effects 
in superconductors were intensively investigated, see for example Ref. ||, §,|5| and 
references therein. 

Multiband superconductors show small values of AC(rc), negative curvature of 
the upper critical magnetic fleld Hc2{T) near the transition temperature, etc. One 
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of the fundamental properties of niultigap superconductors is that nonmagnetic 
impurities are pair breaking ;El for a nice introduction in the properties of muhiband 
superconductors the reader is referred to the review by Moskalenko, Palistrant and 
Vakalyuk.cl Here we consider the case of clean superconductors. 

2. Model and notation 

To begin with, let us introduce standard notations for the dimensionless quasi- 
momentum p and momentum-space averaging in the D-dimensional case, 

^i:/(p)-f'^^f'^/(p) (1) 

where Af is the number of p-points in the momentum-space summation or, equiv- 
alently, the number of unit cells in the crystal. It is also expedient to introduce a 
non-normalized integration over the Fermi surface 

(/(p))F^(5(efc,p-i?F)/(p))p 

where dSb,p is an infinitesimal surface element of the Fermi surface sheet of the bth 
energy band, and v^.p — VpSb^p is the quasiparticle "velocity" which according to 
the present convention has dimension of energy. Conversion to the true velocity in 
m/sec can be performed by multiplying with the lattice constant and dividing by 
h. In this notation for the electronic DOS per unit cell and spin we have 

p(^f) = (1>f- (3) 

We assume that the DOS can be represented as a sum of a regular around the Fermi 
level function pi and a singular part p2 divergent at the energy of the van Hove 
transition -EvhSj 

, , , , , , I pife) = const, £ = i?F, 
p{s)^ Pi{e) + p2ie), i"^)' ' . (4) 

Using the so introduced DOS the well-known formula for the normal specific heat 
per unit cell reads 




cosh2(i/p)^p 

^TT^klT (piiEp) + I P2{Ef + 2k^T v)q,{y)dv ] , (5) 



/ p2{EF + 2k^Tv)q,{v)dv\ 

J — oo / 



where 



£b,p — Ep 

2kBT 



6 / \2 f+°° 
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(see also Fig. |T|). 

3. Ginzburg-Landau coefficients 

The starting point of our thermodynamic analysis is the GL expansion for the 
free energy per unit cell as a function of the temperature T and the order parameter 

— ^=ao^^|-| +2^-1 . (7) 

In thermodynamics of second-order phase transitions the ratio of the GL coefh- 
cients oq and h determines the jump of the specific heat per unit cell at the critical 
temperature, so for superconductors 

where Cg is the specific heat of the superconducting phase. nni— i 
Employing the finite-temperature Bogaliubov-Valatin variational approachQIa' ^ 
we have recently given a simple derivation^ of Gor'kov and Melik-Barkhudarov'stl 
results for the GL coefficients of a clean weak-coupling superconductor with anisotropic 
gap, 

A,,p(r) = S(T)xfc,p. (9) 

The gap anisotropy factor Xb,-p is the eigenfunction of the linearized BCS gap equa- 
tion, corresponding to the maximal eigenvalue. Performing appropriate modifica- 
tion of these results, Eqs. (27, 30) of Rcf. ^ can be written in the form 



|2 



where 



„ , , 1 / tanh V 1 , , , 



Here we have partitioned the bands into a regular and a singular one. For the case 
of regular bands only the general expression for the relative jump of the specific 
heat readsE 



AC 12 1 1 (|ApP)p ^ 12 



7rc^\R ' R /i\ /lA |4\ 3^:7:77 = 1.42613..., (12) 

7C(3)/3a /3a (1)p (|Ap|'i)p 7C(3) 



where the last numerical value is the universal BCS ratio. The question we pose 
now is whether the existence of a VHS can change the inequality (^2|) and thus 
drive an enhancement of the relative specific heat jump, 1//3a > 1- For the latter 
we shall derive a suitable for experimental data processing formula. 
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Figure 1: Plot of the qa (solid line), qh (dashed line), and qc (dash-dotted line) 
functions. 



4. Influence of the van Hove singularity 

We assume that the sample is characterized with sufficient purity, and that the 
VHS is close to the Fermi level, 



n/rc<fc3rc~ I^F-^VHsl 



(13) 



where is the scattering time. In analogy to the normal specific heat (|^), for 
the GL coefficients we can perform the averaging over the Fermi surface for all 
regular bands. For the band having VHS, instead, one has to carry out a separate 
summation over the different constant-energy layers in the momentum space: 



ao = |xiPpi(£'f) + |X2| / p2iEF + 2k^Tc:^)qa{i^)dv, 
u 7C(3) 



/ + 00 
P2{EF + 2k^T,v)qb{v)dv 
-OO 



where 



1X2 I 



lim 

E — >E'p 



\Xi\ 



IX2,p| 



(IXi,pl''%, 



/F,l 



£2,p— 



^^2,p(27r) 



dS- 



2,P 



with fc = 2, 4, and 



2 cosh^ u ' 



S2.p = E 



( tanh I 



W2,p(27r) 



D 



1 



14C(3) V V co^\^v)v-^ 



1 



(14) 



(15) 



7^- (16) 
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The same normalization as for Qc holds for the qa^t functions: j^'^ (la.b{v)dv = 1. 
These functions are shown in Fig. If, in an acceptable approximation, we can 
consider the gaps in the two different bands to be nearly constant, the sign of the 
energy-surface averaging can be dropped: 



P2{ET. + 2k^T,v)qa{v)dv, 

-OO 

^^^^^ ^|Xi|Vi(i?F) + |X2|' P2{E^ + 2KT,v)q,{v)dv 

J — OO 



(17) 



The generalization of these expressions for the multiband case is straightforward. 
Using (ITh, we find for the specific heat jump 



AC 

cT 



12 1 



7C(3) /3vHS ' 



(18) 



where for the renormalizing multiplier we obtain a generalized two-band expression, 
cf. Refs.|,|ll,| 



/ + 00 
-OO 



P2{.£iy)qc{v)dv 



/ + 00 
-OO 



v)dv 



with 



(19) 



(20) 



In the multiband case a summation over the band index is required, and the gener- 
alized Moskalenko-Pokrovsky formula reads 



AC 

C„ 



12 
7C(3) 



Pb{eu)qa{v)dv 

/ + 00 
Pb{eu)qc{v)dv 
-OO 

/+00 
Pb{e^)qb{v)dv 
-OO 



(21) 



Note that the GL order parameter ^ in (|l9|) cancels. This is related to the condi- 
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Figure 2: Typical constant energy surfaces for cubic perovskites.E^ When the "dices" 
(left) become cubes the perfect nesting of their flat surfaces creates ID singularity 
(oc Xj \Je — i?vHs) of the DOS. Similarly, the nearly constant cross-section of the 
narrow tubes (right) creates a 2D singularity of the DOS. The box indicates the 
first Brillouin zone. 



tional factorization 

= ^Xp = (cS)(xp/c), 
aoiSp = (ao/c2)|cSp « IxpPlSp = {\xp\^ j ^)\cE\\ 
= (Vc')|cS|4 (X IxpllSl^ = (|xp|Vc^)|cS|^ 
AC _ \_al _ 1 (flo/c^)' 

^ |Sp + ifel-r = (ao/c^)^ IcSp + \{h/c%c-\\ (22) 

which conserves as scaling invariants all measurable quantities like, e.g., the heat 
capacity and the free energy. Using this freedom we can normalize |xpP — 1- 

Analyzing (|l^) one can easily verify that for the model cases of one-dimensional 
(ID) and two-dimensional (2D) VHS, respectively, 

p(^°)(e) = (e - i^vHS )-'/', £ - i^vHS > 

p('°)(£) = -ln|£-i?vHs|, (23) 

the renormalizing multiplier /3yyg could be > 1. In this case the VHS can enhance 
the relative jump AC/Cn|Tc to values bigger that the conventional 1.43. For getting 
an insight into the influence of the VHS we recommend the DOS p2{Ep + 2k^Tc v) to 
be plotted along with the qa,b,c{v) functions. There is no mathematical restriction 
how big l//3vHS could be for a very narrow peek of the DOS. However, realistic 
DOS corresDonds to ID or 2D case. Such type of singularities appear for cubic 
perovskitealj where close to the VHS some constant-energy surfaces have the shape 
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of elongated tubes (2D cross section), or dices (ID cross section). Some typical 
surfaces are shown in Fig. |^. Just the opposite situation could arise, however, in a 
two-band VHS model if the Ep lies in between two logarithmic VHS. In this case Cn 
can be considerably enhanced. We find it instructive for a realistic VHS model the 
thermodynamic variables, e.g., AC/Cn(Tc), AC and Cn(rc), to be plotted versus 
E-p — Eyhs- The Cn{T)/T dependence for some particular {Ep — £'vHs)-values is of 
general interest even far from T^. 

Concluding, we believe that the influence of the VHS can be studied by fo- 
cusing on a very simple characteristic — the jump of the specific heat at Tc- Such 
measurements do not require large single crystals. 
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